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Abstract 

We present a detailed comparison of two self-consistent equation-of- 
state models which differ from their electronic contribution: the atom 
in a spherical cell and the atom in a jellium of charges. It is shown 
that both models are well suited for the calculation of Hugoniot shock 
adiabats in the high pressure range (1 Mbar-10 Gbar), and that the 
atom-in-a-jellium model provides a better treatment of pressure ion- 
ization. Comparisons with experimental data are also presented. Shell 
effects on shock adiabats are reviewed in the light of these models. 
They lead to additional features not only in the variations of pressure 
versus density, but also in the variations of shock velocity versus par- 
ticle velocity. Moreover, such effects are found to be responsible for 
enhancement of the electronic specific heat. 

1 Introduction 

The studies concerning matter under extreme conditions have broad ap- 
plications to material science, inertial confinement fusion and astrophysics. 
About 98 % of the Universe is made of hot and ionized matter. The center 
of stars and white dwarfs have pressures exceeding thousands of Mbar and 
temperatures of tens of millions of K. In neutron stars pressures of millions 
of Mbar are believed to exist. 

Nowadays, the shock-wave technique in laboratory laser experiments [1] 
enables one to reach maximum pressures that are as high as hundreds of 
megabars and even more for many materials |2]. Such pressures range to 
three orders of magnitude higher than the pressure in the Earth center and 
close to that in the center of the Sun. 

Therefore, the need for suitable equation of state (EOS) of high-energy- 
density matter becomes crucial. The thermodynamics and the hydrodynam- 
ics of these systems can not be predicted without a knowledge of the EOS 
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which describes how a material reacts to pressure. For instance, the theory 
of stellar evolution is affected by uncertainties in EOS (the internal structure 
and cooling time of the white dwarfs depend on the details of their EOS). 
More precisely, in astrophysical objects such as low-mass stars, brown dwarfs 
and giant planets, hydrogen and helium highly dominate, and their EOS, 
including non-ideality effects, can be accurately evaluated using "chemical" 
free-energy models (see for instance 0H1E]). However, the EOS of carbon 
is also required for modeling inner envelopes of carbon-rich white dwarfs 
or outer accreted envelopes of neutron stars [Hj. For general astrophysical 
applications of various types of stars, the EOS of oxygen is also required. 
In the center of the Sun, where matter density is about 157 g/cm 3 and 
temperature 1.3 keV, although helium and hydrogen are the most abundant 
elements, the small quantity of iron has a strong impact on the calculation 
of the radiative opacity. Therefore, it is important to be able to determine 
the partial density of iron, which requires knowledge of its EOS. 

In the present work, we consider strongly coupled (nonideal) plasmas, 
characterized by a high density and/or a low temperature. In such plas- 
mas, ions are strongly correlated, electrons are partially degenerate, the De 
Broglie wavelength of the electron becomes of the same order of the inter- 
particle distance, and the coupling parameter 



ratio of Coulomb potential energy and thermal energy, is much larger 
than unity. The matter density re, represents the number of particles per 
volume unit, Z* the average ionization of the plasma, ks the Boltzmann 
constant and T the temperature (throughout this paper, physical quantities 
are expressed in atomic units (a.u.)). 

When a material is subjected to a strong shock wave, it becomes com- 
pressed, heated and ionized. As the strength of the shock is varied for a 
fixed initial state, the pressure-density final states of the material behind 
the shock belongs to a curve named shock adiabat or Hugoniot curve. The 
Hugoniot curve depends on the EOS of the matter, which, in principle, can 
be determined from theory. In practice, the barriers to ab initio calculations 
are formidable owing to the computational difficulty of solving the many- 
body problem. Consequently, it has proven necessary to introduce simplify- 
ing approximations into the governing equations. The price of simplification, 
however, is typically a loss of generality and the resulting theory is adequate 
only over a limited region of the phase diagram. The determination of the 
average electronic charge density relies usually on the Density Functional 
Theory. A well-known example is the Thomas-Fermi (TF) model of dense 
matter [7], which contains certain essential features in order to characterize 
the material properties over a limited range of conditions. At intermediate 
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shock pressures, when the material becomes partially ionized, the EOS de- 
pends on the precise quantum-mechanical state of the matter, i.e. on the 
electronic shell structure [BJ[8]. Therefore, there is a great interest in the 
physics of bound levels in high-energy-density plasmas [9] and quantum self- 
consistent-field (QSCF) models are replacing the TF approach. Such statis- 
tical models are well-suited for elements with a few occupied orbitals, i.e. 
with a "reasonable" number of electrons. Therefore, the theory presented 
in this paper is not applicable to hydrogen and helium, since their atomic 
number is too low. A lot of theoretical and computational efforts have been 
made in the United States and in Russia in the last twenty years concerning 
quantum effects on Hugoniot shock adiabats. It is important to review that 
huge amount of work, in the light of the most recent developments in the 
modeling of atomic structure. Indeed, QSCF calculations strongly depend 
on the modeling of the ions and on the treatment of free electrons. We 
present a QSCF-based EOS model, which will be named ESODE (Equation 
of State with Orbital Description of Electrons) in the following. The ionic 
contribution to the EOS is described by a perfect ideal gas, and the cold 
curve (T=0 K isotherm) has been obtained in most of the cases from Aug- 
mented Plane Wave (APW) calculations [TO]. Exchange-correlation effects 
at finite temperature are taken into account. ESODE comes in two main 
versions for the thermal electronic contribution to the EOS: Average atom 
in a Spherical Cell (ASC) or Average atom in a Jellium of Charges (AJC). 
In both cases, bound electrons are treated quantum-mechanically. In the 
ASC model, all electrons are confined within a Wigner-Seitz sphere, and 
the continuity theorem of the wavefunctions at the boundary can be satis- 
fied in different ways. In the AJC model, bound-electron wavefunctions can 
extend outside the sphere, where the plasma is represented by a uniform 
electronic density (jellium) neutralized by a continuous background of pos- 
itive charges, representing ions. An interesting feature of AJC model is the 
possibility to use a quantum-mechanical description of free states (instead 
of a semi-classical one) together with a careful search for shape resonances. 
This seems to be necessary in order to ensure thermodynamic consistency, 
and a suitable description of polarization and density effects. The QSCF 
models (ASC and AJC) are described in section[2l Shock adiabats calculated 
from those models are presented, analyzed and compared to traditional TF 
model and to published experimental data in section El An interpretation 
of shell effect in terms of specific heat is discussed in section [3] as well. The 
dependence of shock velocity on particle velocity is analyzed in section [H 
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Figure 1: Hugoniot curves for Be. prj =1-83 g/cm 3 . Experimental values are 
given by Refs. [33 M M ESI HE EZ] • 

2 The models 

2.1 First quantum statistical model: Atom in a Spherical 
Cell (ASC) 

The task is to evaluate the contribution to the EOS coming from the exci- 
tations of the electrons due to temperature and compression. Atoms in a 
plasma can be idealized by an average atom confined in a Wigner-Seitz (WS) 
sphere, which radius r ws is related to matter density. Inside the sphere, the 
electron density has the following form 



n(r) 



7T" 



+X^3/2 {-V(r),fi-V(r),x )], 



(2) 



where 



fi{x,y) 



2(21 + 1) 



1 + exp[(x - y)/k B T] 
is the usual Fermi-Dirac population and 



(3) 



4 



~ [°° n/2 (1+ /2) l/2 

J n/2 (a,x,a) = 1 + ej/ _ s (4) 

is the modified Fermi function of order n/2. The first term in corre- 
sponds to the contribution of bound electrons to the charge density, while 
the second term is the free-electron contribution, written in its semi-classical 
TF form. The energy e;, and wavefunction ip b of a bound orbital are calcu- 
lated in the Pauli approximation j!2j . in which only first-order relativistic 
corrections to the Schrodinger equation 

- ^AV>6 + [V(r) + V mv (r) + V D (r)]^ b = e b ^ b (5) 

are retained. We introduce the notations % = ^bT/Eq, Eq being the 
rest mass energy of the electron and V(r) = V(r)/(k B T), where V(r) is the 
self-consistent potential: 

V(r) = + / -^-UV + V xc (r), (6) 

r Jo \r — r'\ 

V xc being the exchange-correlation contribution, evaluated in the local 
density approximation [13]. V mv is the mass- velocity correction 

and Vd the Darwin correction 

I? / \ 1 1 dVfd In 

V ° {r) = -2 2E + e b -V(r) ^Kd-r ~ rY (§) 
Omission of the spin-orbit term is a consequence of the spherical sym- 
metry. Last, the chemical potential fi is obtained from the neutrality of the 
ion sphere: 



pi 
Jo 



n(r)47rr dr = Z, (9) 

10 

and p = n/(k B T). Eqs. © and © must be solved self- 

consistently provided that, at each step of the iterative process towards 
convergence, bound orbitals are obtained from the Schrodinger equation. 
The electronic pressure P e [HI [15] consists of three contributions, P e = 
P b + Pf + P xc , where the bound-electron pressure P b is evaluated using the 
stress-tensor formula 



+ 2e 6 (l + 
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Figure 2: Hugoniot curves for B. po =2.34 g/cm 3 . Experimental values are 
given by Ref. [321 ES]. 

Ub representing the radial part of the wavefunction ipf, multiplied by r. 
The free-electron pressure Pf reads 

Pf = ^(k B T) 5 / 2 [J 3/2 (-V(r ws ),fi-V(r ws ), X ) 

+ ^J 5 /2(-V(r ws ),fl-V(r ws ),x)] (11) 

and P xc is the exchange-correlation pressure evaluated in the local den- 
sity approximation [13]. The choice of the boundary conditions plays a major 
role in the expression of pressure. We chose a decreasing-exponential bound- 
ary condition, since it allows the matching of wavefunctions with Bessel 
functions outside the WS sphere [Hj, where the potential is zero. Moreover, 
such a condition is consistent with the fact that the bound-electron density 
vanishes at infinity. The internal energy in the ASC model is 



E e = y^njej-- r W \(r) [ TWS ^Xd 3 rd 3 r' 

rr ws 

+E XC - / n(r)V xc {n{r))d 3 r, (12) 



where E xc is the exchange-correlation internal energy and tia the popula- 
tion of state j (either bound or free). The first term in (I12p can be expressed 
by 
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Figure 3: Hugoniot curves for C (graphite), po =2.28 g/cm 3 . Experimental 
values are given by Ref. [33 E3 ESI S3] ■ 



(13) 



Ef tP is the potential energy 



E 



f,P 
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[Ji/ 2 (-nr),M-F(r), X ) 



+xJ 3 /2(-V(r),fl - V(r), xW(r)d 3 r 
and Ef t k the kinetic energy 



(14) 



E 



f,k 



V2(k B T) 5 / 2 



7T- 



[^3/2(-V-(r),A-^(r),x) 



+xJs/2(-V(r),fi-V(r),x)]d 3 r. 



(15) 



In the present work, only non-relativistic calculations are carried out, 
which correspond to Eq — ► oo in Eqs. ([2]), (flOj) . (fTTj) . (fT4"l) and (fl~5l) . in order 
to make comparisons with the model described below (in section [2~2"j) . which 
is non-relativistic. Formulas ([2]) to ()15p enable one to include relativistic 
effects without solving Dirac equation, in which the wavefunction consists 
of two components, making the choice of boundary conditions even more 
awkward. 
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Figure 4: Hugoniot curves for Al. po =2.70 g/cm 3 . Experimental values are 
given by Ref. [33 SHI SH S3 H3 ISl S3 SHI [13 SS] - 

2.2 Second quantum statistical model: Atom in a Jellium of 
Charges (AJC) 

In order to go beyond the TF approximate treatment of continuum electron 
charge density, it is necessary to use a full quantum-mechanical description 
of the continuum states and to consider that both bound and free orbitals 
can extend outside the WS sphere. This leads us to define the environment 
outside and far from the central average ion. One way to address this is the 
jellium model (or electron gas model), i.e. a uniform electron density — n + 
neutralized by a positive background n + simulating the ionic charges. The 
model relies on a method proposed by J. Friedel [TBJUZ] in order to treat the 
electronic structure of an impurity, represented by a spherical potential of 
finite range, in an electron gas. It has been further developed by L. Dagens 
[181 [T9] and F. Perrot [201 |2T| 122j . In this framework, the electron density 
can be written: 

n(r) = -n+ + An a + An ab + An abc + (16) 

where — n + is the density of the electron gas (jellium), An a the charge 
density displaced by the immersion of an atom "a" in the jellium, An a h the 
difference between the charge density of the a — b molecule and the charge 
densities of free atoms a and b, etc. The AJC approach retains only the 
first two terms. A nucleus of charge Z is then introduced in a cavity in the 
positive background. The radius of the cavity represents the average atomic 
radius in the plasma. Therefore, the problem is reduced to the response 
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Figure 5: Hugoniot curves for Fe. po =7.85 g/cm 3 . Experimental values are 

given by Ref. 0aiSBISlHai5DlEIll521l53ll53ll551l5aEZI- 



of the electrons to the immersion of a positive charge Z. Ionic response is 
roughly simulated by the formation of a cavity. These considerations lead 
to the following form of the electron density: 



2 r°° 



de, 



where the eigenstates (ipb> e b) an d the scattering states (tft e j,e) are ob- 
tained solving Schrodinger equation with the potential 



V(r) 



Z 

r ./u 



H 



n(r f ) 



d 6 r' + V xc (r)-V xc (-n+), 



(17) 



where Roo » r a ,r a being the radius of the cavity. As in the ASC model, 



V xc is calculated from [13]. V(r) is determined, as in the ASC model, in a 
self-consistent way. The ionic density is modeled by 



n+(r) = 

The expression of internal energy is [22 



for r < r a 
n + for r > r a . 



(18) 



E[n,n+] = Z*[e k + e xc ](-n+) + AE[n,n + ] - A T Piv a , (19) 
with v a = 4vrr 3 /3, A T = (dlnZ* /d\nT)\ VajT and 
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Figure 6: Hugoniot curves for Cu. po =8.84 g/cm 3 . Experimental values 
are given by Ref. ^EilSB|3a|3giSll5Sll521ISaEIlE2lE31IS3lE51. 



Pi = — f ^r 2 V(r)dr. (20) 

V a Jr a 

e k and e xc are the kinetic and exchange-correlation energies per free 
electron and AP[n,n+] is the energy change resulting from the immersion 
of an ion in the jellium: 

f Ro ° - 3 

AE[n,n + ] = / (e k [n(r)) - e k [-n + ])d r 

Jo 

Hoc 2? 

— (n(r) + n + (r)W 3 r 
r 

2 Jo Jo \r — r'\ 

[■Roc, 

+ / (e xc (n(r)) - e xc (-n+))d 3 r. (21) 
Jo 

The pressure is obtained by 

P e = [P k + Pxc](-n+) + n+^(r a ) + (1 - A^,t)Pi, (22) 
where A„ a = (dlnZ* /dlnv)\ Vaj T and 
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Figure 7: Internal energy versus compression along the Hugoniot curve for 
Be, B, Al and Fe in the ASC model. 



The term Pk corresponds to the pressure of a free-electron gas. The 
expression of pressure in (|22p is rigorously the derivative of energy versus 
volume. This is the main difference with INFERNO model |23} 124] . The 
major difficulty of these models is that the average ionization is not well 
defined when the outer electrons are more or less delocalized. The only 
way to make the formalism variational is to specify the ionization in the 
AJC model. In other words, the question is how to define the residual 
electron density (— n+) far away from the point where the positive charge 
is introduced into the jellium. A convenient choice is Z*(v) = Z^ F (v) |22j . 
Z?pp being the TF ionization. In such a way derivatives of ionization with 
respect to volume and temperature can be obtained analytically, using the 
numerical fit proposed by R.M. More |25} 126]. 

ASC and AJC models are very different concerning the modeling of the 
environment of the atom, isolated and confined in the ASC model, and 
immersed in an infinite effective medium in the AJC model. 




(23) 
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Figure 8: Compression versus temperature (thermodynamic path) along the 
Hugoniot curve for Be, B, Al and Fe in the ASC model. 

2.3 Nuclear contribution and cold curve 

The adiabatic approximation is used to separate the thermodynamic func- 
tions into electron and nuclear (ionic) components. The total pressure can 
be written: 

P(p,T) = P c (p) + Pi(p,T) + P t (p,T), (24) 

where Pt(p, T) = P e (p, T)—P e (p, 0) is the thermal electronic contribution 
to the EOS, P e being electronic pressure. Pi is the nuclear (ionic) pressure 
and P c the cold curve. The replacement of P e (p, 0) by P c (p) is necessary 
since the results of the QSCF models are not valid for T — > 0. Following the 
same procedure, internal energy is divided into electronic, ionic and cold 
component. Expression (|24p also holds for internal energy. The nuclear 
component is evaluated in the ideal-gas approximation. 

• The cold curve must satisfy the condition that the quantity 

p{dP/dp) (25) 

evaluated at normal density p n must be equal to the experimental bulk 
modulus Bq, which can be obtained from the sound velocity c s through 
B Q = pc 2 s . The cold curves are obtained either from APW [TU] simulations, 
or using the Vinet [27] universal EOS. In the latter case, pressure can be 
evaluated analytically by 
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Figure 9: Absolute value of eigen-energy of orbitals Is, 2s and 2p of Fe 
versus compression along the Hugoniot in the ASC model. For a given value 
of p, the corresponding value of T can be found on Fig. El 



P c = [35 (l-x)/x 2 ]exp[C(l-x)] (26) 
and the internal energy by 

E < = - ex PtC(l - *)][! " C(l - *)]] (27) 

with C = 3(B' Q - l)/2, x = (p/p n ) 1/3 and B' = (dB /dP)\ P= o,T- Values 
of such parameters obtained from experiments are available for several ele- 
ments in the literature (see for instance [27\ I28j). In many situations, Vinet 
EOS gives realistic results, and its accuracy is sufficient for most applications 
involving high pressures and temperatures. 

3 Shock Hugoniots 
3.1 Definitions 

The initial state of the plasma is characterized by a density po, a temperature 
To, a pressure Po, and an internal energy Eq. D is the shock velocity, u the 
matter velocity, P and E are respectively the pressure and internal energy 
behind the shock front. Since the mass preservation has to be satisfied, the 
density of the compressed gas p satisfies the relation poD = p(D — u). The 
resultant force acting on the compressed gas is equal to the difference of the 
pressure on the shock side and on the side of the undisturbed fluid, that is 
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Figure 10: Electronic specific heat for Al. 

P — Pq = pqDu. The increase in the sum of the internal and kinetic energies 
of the compressed gas is equal to the work done by the external force acting 
on the shock front, i.e. pqD(E — Eq + u 2 /2) = Pu. Rearranging these 
equations, one gets [29] : 

hp + P )(±-l)=E-E , (28) 

which is known as the Hugoniot relation. It is worth keeping in mind that 
the Hugoniot shock adiabat is not the collection of successive states of the 
matter during the propagation of the shock. The only relevant information 
is that the state of the matter "after the shock" belongs to the Hugoniot. 
Technically, the Hugoniot curve can be obtained by solving (|28p at each tem- 
perature step, either calculating the EOS "in line" or by an interpolation in 
a table (pre-tabulated EOS). The two relations P(p,T) and E(p,T) consti- 
tute the EOS, which thermal part is calculated from the models described 
in sections 12.11 and 12.21 

In this work, focus is put on the standard Hugoniot (Po=0, po equal to 
solid density and To=300 K). No distinction between po and p n will be made 
in the following. As mentioned before, total energy is preserved during the 
shock, since an increase of internal energy is transformed into heat, which 
makes the temperature of the plasma increase. More exactly there is no heat 
exchange during the shock (which explains the name "shock adiabat"), but 
an increase of internal energy, distributed in internal degrees of freedom. 
One can wonder whether this conservation of energy is in contradiction 
with the fact that a physical system tends to minimize energy, as invoked in 
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Figure 11: Electronic specific heat for C. 



sectionEJ In classical mechanics, Noether's theorem, which links symmetries 
and conservation laws, relates the invariance by translation in time to the 
conservation of total energy. Therefore, the idea that a physical system tends 
to minimize energy is often included in models where time is absent. This is 
also the case of the present study, since in QSCF models the wavefunctions 
are assumed to be time-independent. Therefore, only stationary solutions 
are considered here, and only systems in thermal equilibrium are supposed to 
exist. There is actually no contradiction: total energy is preserved during the 
shock, but the pressure calculated in the AJC g model is thermodynamically 
consistent since it is calculated as a derivative of the energy (see (I22p in 
section [272]) . Calculations have been performed for T < 6.4 keV. Figures [TJ 
El El HIE] and [6] represent Hugoniot curves for beryllium (Be, Z=4), boron (B, 
Z=5), carbon (C, Z=6), aluminum (Al, Z=13), iron (Fe, Z=26) and copper 
(Cu, Z=29) respectively, calculated from pure Thomas-Fermi EOS, ASC 
model, and AJC model with a full quantum treatment of electrons (namely 
AJCq). The calculation within the jellium model with a hybrid description 
of electrons, i.e. where bound electrons are described in the framework 
of quantum mechanics and free electrons in the TF approximation, is also 
presented (AJC/i), since it is not equivalent to ASC model, even if it relies 
also on a hybrid description of electrons. The cold curve has been calculated 
using APW method for Be, C, Al, Fe and Cu. For Be, B, Al and Fe, 
Vinet universal EOS gives, for the Hugoniot curve, results which we find 
sufficiently accurate. For B, the Vinet parameters have been interpolated 
between Li (Z=3) and Be. All the QSCF models give very close results 
and differ from the TF approximation; it seems however that both hybrid 
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Figure 12: Electronic specific heat for Fe. 

models (atom in a spherical cell ASC or in a jellium of charges with quantum 
bound electrons and TF free electrons AJC/J are very close to each other 
and differ from the full quantum atom-in-a-jellium-of-charges model AJC g . 
For C, the Hugoniot curve obtained from ASC model is different from the 
ones obtained from other QSCF models, since the K shell is ionized for a 
lower density in the ASC approach, which induces a larger pressure. In the 
case of Al, the difference between the theories appears for P > 3 Mbar and 
AJCg gives higher pressures than the other models for 2< rj < 3.5, where 
r\ = p/po is the compression rate. For r] > 3.5, all QSCF models give lower 
pressures than TF model. 

3.2 Quantum shell effects 

Our models emphasize the thermodynamic domain where the Hugoniot 
curve strongly depends on the electronic structure, i.e. beyond four times 
solid density where the shoulders (double in the case of Al, Fig. 2|) corre- 
spond to ionization of successive shells in the Average-Atom picture. The 
Hugoniot curve tends to the classical limit equal to four times the solid 
density at very high pressure, or typically T > 10 keV. These shoulders 
are a consequence of the competition between the release of energy stocked 
as internal energy within the shells and the free-electron pressure. When 
ionization begins, the energy of the shock is used mainly to depopulate the 
relevant shells and the material is very compressive. However, the pressure 
of free electrons in increasing number dominates again and the material be- 
comes more difficult to compress. Both models show compression maxima 
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Figure 13: Shock velocity D versus particle velocity u for Be, C, Al and Fe. 

in the range 5po — 6po. In this region, the electrons from the ionic cores are 
being ionized and the shock density increases beyond the infinite pressure 
of 4po- As ionization is completed, the plasma approaches an ideal gas of 
nuclei and electrons and the density approaches the fourfold density 4po . 

The successive electron shells are accurately represented. For Be (Fig. 
[1]), B (Fig. [2]) and C (Fig. [3]), all the models show a single density maximum, 
corresponding to the ionization of the K electron shell. For Al (Fig. H|), 
there are two density maxima corresponding to the K and L electron shells. 
For Fe (Fig. [5]) and Cu (Fig. [6]) there are density maxima or inflexions 
corresponding to the K, L and M electron shells. The L shell ionization 
feature gives the largest density increase. In the case of Cu (Fig. [6]), the 
ASC and AJC^ models exhibit a kind of discontinuity around 40 g/cm 3 , due 
to pressure ionization of 2p orbital. Such a sharp increase of pressure does 
not exist in the AJC g model, because of its thermodynamic consistency. 
This is due to the fact that it relies on a variational formulation: pressure is 
rigorously obtained as a derivative of the free energy, and shape resonances 
are carefully taken into account in the quantum treatment of free electrons. 
Such features lead to a continuous disappearing of a bound state into the 
continuum [30]. This is not the case of ASC and AJC/j models. In the 
framework of these two models, we tried to improve the dissolution of bound 
states into the continuum splitting them into bands and using a simple 
degeneracy reduction; this makes the discontinuity a little less sharp but 
does not solve the problem. Only a careful search of resonances would 
help, together with a variational calculation of pressure. Figure [7] illustrates 
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29 Cu 


6.846 


5.905 


1.137 


1.212 


6.88 10~ 4 


40 Zr 


6.926 


/ 


1.125 


/ 


3.28 10" 4 


42 Mo 


6.975 


6.711 


1.123 


1.149 


3.50 10" 4 


4 8 Cd 


3.676 


4.251 


1.180 


1.182 


9.07 10~ 5 



Table 1: Shock adiabat parameters of Be, B, C, Na, Mg, Al, Fe, Co, Cu, Zr, 
Mo and Cd. Coefficients 00 and a± are compared with those published for 
metals in Ref. |67] . 



the fact that internal energy also reflects the oscillations due to the shell 
structure, as well as the thermodynamic path (T, p) (see the cases of Be, B, 
Al and Fe in Fig. [8]). In fact, the non-monotonic character of thermodynamic 
variables stems from the eigen-energies of the orbitals themselves, which 
exhibit the oscillations as well (see Fig. [9]). The first density for which a 
lower compression is obtained is named "turnaround" point. 

The pressure differences from the quantum mechanical theory along with 
the oscillatory feature at the turnaround density, can be explained by ex- 
amining the heat capacity predicted by the TF and ASC theories along the 
Hugoniot path. At low temperature, the electronic heat capacity depends 
on the number of electrons that can be excited around the Fermi energy. 
The TF theory predicts a smooth increase since the density of states in this 
model is a monotonic function of energy. Therefore, the electronic specific 
heat 

r elec _ _J_ d[E(p,T) - Ej(p,T)] 

v ~ 3k B /2 8T { yj 

is an interesting indicator of the ionization of successive shells and indi- 
rectly on the release of energy. Figure [TO] represents Cf ec along the Hugo- 
niot curve in the TF and in the ASC models for Al. Both theories show the 
effect of the coulomb attractive potential of the nucleus binding the elec- 
trons, represented by the peak around 300 eV for TF theory and 100 eV for 
ASC model. After the turnaround, there are 11 free electrons and 2 bound 
electrons remaining in the Is orbital (K shell), which is far away from the 
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Figure 14: Derivative of shock velocity dD/du versus particle velocity u for 
Be, C, Al and Fe. 

energy zero (1.5 keV at T=100 eV). This phenomenon is a kind of Schottky 
anomaly [31J. As long as temperature is not sufficient to ionize those two 
electrons, the specific heat tends to an asymptote corresponding to an ideal 
gas of 11 independent particles. When both Is electrons are ionized (after a 
"threshold" temperature), there is a sudden break in the specific heat, which 
tends to an ideal gas of 13 independent particles. It is not as important for 
the 2s and 2p bound states (L shell), since their energy levels are not as far 
from the continuum (a few tenth of eV). Figure [TT1 represents the electronic 
specific heat along the Hugoniot curve in the TF and in the ASC models 
for C. In that case, the only difference between both models is that there is 
a strong enhancement due to ionization of Is orbital in the ASC approach, 
followed by the asymptotic limit of 6. The same phenomenon occurs with 
Fe (Fig. [T2|) ; in that case, after two Schottky anomalies, the electronic part 
of the specific heat tends to an ideal gas of 26 electrons. 

3.3 Comparisons with experimental data 

Experimental data on Be [331 El US [361 EZ], B [35], C [351 EH EH], Al 

[iQl SH S21 SH HH SSI HH S71 HE] , Fe ^HISII3SII5QIEHI521E31E31I551ESIE23 

aiidCu|3aiiDll^li3l^|^l5aiSIlEIlE2lE3llS3lE51 have been collected 
for comparison with the EOS models presented above. Several experimen- 
tal methods have been used to generate well-defined shock states: gas guns 
for pressures up to 5 Mbar, explosive-driven spherical implosions and laser- 
driven plane waves for generating shocks up to 10 Mbar. High-power lasers 
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are capable of driving multimegabar shocks in small samples either by di- 
rect irradiation or indirectly. In the latter case, the target is mounted on 
a hohlraum into which the lasers are focused. Pressures up to a few Gbar 
have been reached through underground nuclear explosions. If the shocks 
have been measured accurately, superimposing all of these data should yield 
a single smooth shock Hugoniot curve. The maximum pressures reached in 
the experiments are: 18 Mbar for Be, 4000 Mbar for Al, 191 Mbar for Fe 
and 204 Mbar for Cu. It appears difficult to decide which model gives the 
best agreement with experiments, since there are obviously very few avail- 
able data for the region of interest (typically above 100 Mbar). The error 
bars associated to the ultrahigh-pressure experimental values for Al are too 
large to conclude concerning the existence of shell effects. For Fe and Cu, 
results from AJC q model seem to agree better with the first interpretable 
experimental points than the hybrid models (ASC or AJC/J. The large er- 
rors in the measurements make the experimental results useless for refining 
the theoretical data. Therefore, it is not really legitimate to attempt to re- 
late gas-dynamics measurements with a discussion of shell effects. Anyway, 
analysis of the computational results shows that the deviation from the ex- 
perimental points of Al can not be explained only by shell effects. In the 
region where most of the experimental points are available, the role of the 
cold component is important, while shell effects begin to play a significant 
role after the matter is already compressed (around the asymptotic com- 
pression 77 ~ 4) and begins to heat up. In order to discriminate models from 
eachother, the uncertainty on velocities should be very small (typically < 
1%), which is not possible. However, if the temperature beyond the shock 
could be measured with an incertainty of about 20 %, it should be possible 
to decide which model gives the best agreement with experimental data in 
the (T, P) plane. 



4 Shock and particle velocities 

The particle and shock velocities read (see section 13. ip : 



it 



' (p-poXP-Po) 
Pop 



and D 



' p(p-Pq) 

Pa(p-po)' 



(30) 



Expressions in (I30p are generic; they do not involve, a priori, any explicit 
relation D(u). However, one has [66J: 



An 



u = Co H u H 

2 6c 



An rG. 



1 P 



(3) 



POP 



(2) 



4 A ° 



u 2 + 0(u 3 ), 



(31) 



where Go is the Griineisen coefficient at p = po, 
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Figure 15: Second derivative of shock velocity d 2 D/dv? versus particle ve- 
locity u for Be, C, Al and Fe. 



A = -P^/(2 P0 P^) (32) 

and 

p (n) = (dWp/d{l/p)W)\ s . (33) 

For a perfect ideal gas, the Griineisen parameter (dP/dE)\ p / p tends to 
2/3. Expression (f3Tj) is suitable for densities close to the pole (po, To, Pq, Eq), 
since it comes from a Taylor development of D and u around that specific 
initial point. For metals, it is often invoked that the relation is quasi-linear. 
Figure [p3] illustrates the fact that the (u, D) relationship is almost linear 
over a wide range of densities, except close to the origin, where we could 
not perform the calculation. It is interesting to try to fit the (u, D) relation 
with a quadratic expression [67] D = a2U 2 + a\u + ao- The parameters are 
presented in tabled] for Be, B, C, Na, Mg, Al, Fe, Co, Cu, Zr, Mo and Cd, 
and compared to the values published by Al'tshuler et al [67J. The slope is 
different from 4/3, corresponding to the perfect ideal gas, as can be checked 
from (I30p . However, if one looks more carefully at the first and second 
derivatives (Figs. Q3] and [To"]) of shock velocity versus particle velocity, one 
finds that the behaviour of shock velocity is more complicated, and that 
there are some oscillations, reflecting the shell structure as well, and even 
inflexion points (see Fig. [To]) . The amplitude of oscillations in the (u,D) 
relationship is very small, which is not the case of (p, P) relationship. This 
is even more sensitive at the turnaround of the shock adiabat. Indeed, in 
that case, (dP/dp)\s has strong variations, which implies, according to the 
relation 
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Figure 16: (D — u) versus u for Be, C, Al and Fe. 



P 



,dP_ 

dp 



+ P^ 

S dE 



9 



,dP_ 
dp 



(34) 



that P has strong variations as well, since (dP/8E)\„ and {dP/dp)\E 
have smooth variations. A survey of Figs. [T3"l [TH and [T^] shows download 
curvature, upward curvature and more complicated deviations from non- 
linear behaviour. G.I. Kerley [US] suggests that it is more convenient to 
plot {D — u) versus u instead of D versus u. His analysis concerns mostly 
low velocities, typically below 10 km/s, but Fig. [TBI shows that it remains 
relevant for high velocities in the region affected by shell effects and confirms 
our results. 



5 Conclusion 

The application of shock waves to plasma physics makes possible the gener- 
ation under laboratory conditions of extremely high energy densities typical 
of astrophysical objects such as stars and giant planets. The physical infor- 
mation obtained this way extends our basic knowledge on physical properties 
of plasma to a broad area of the phase diagram up to pressures nine orders 
of magnitude higher than the atmospheric one. 

We presented and compared two equation-of-state models in which the 
thermal electronic contribution is evaluated using two different quantum self- 
consistent-field models: the atom in a spherical cell (ASC) and the atom in a 
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jellium of charges (AJC). It was shown that both models are very well suited 
for the calculation of high-pressure Hugoniot shock adiabats and that the 
AJC model, relying on a full quantum treatment of the electrons, provides 
a better treatment of pressure ionization. 

We proposed a qualitative and quantitative study of quantum orbital 
effects on the principal shock adiabat for different elements, In the (p,P), 
(/?, E) and (T, p) representations, such effects lead to oscillations correspond- 
ing to the ionization of successive orbitals. We found that such oscillations 
are also visible in the electronic specific heat, through Schottky anomalies, 
and in the energies of the orbitals themselves. They are responsible for devi- 
ations from linear behaviour in the relationship between shock velocity and 
particle velocity. 

The next step will be to check whether the oscillations still exist "be- 
yond" the Average Atom model, i.e. if real electronic configurations are 
taken into account |14t 115] . The remaining difficulty of our models comes 
from the impossibility to represent in a simple and suitable way the influence 
of the plasma environment on a specific ion. Indeed, such an environment 
fluctuates; the number, the localization in space and the structure of neigh- 
bouring ions may change drastically. In principle, it should be necessary to 
calculate a large number of geometric configurations of the system (positions 
of ions in space), and then their thermodynamic average. Such an approach 
is untractable ab initio, except for very low-Z atoms in a restricted range of 
temperature and density, using molecular-dynamics simulations [69j. 
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